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In this paper, we present results of a feasibility study of a model-independent way to measure 
the angle 03 of the unitarity triangle. The method involves — > DK^ decays where the neutral 
D decays to the KgTT^Tv~ final state, together with the sample of decays of CP-tagged D mesons 
(produced, e.g. in ip{3770) — » DD process) to the same final state. We consider different approaches 
to the extraction of (j>3 and obtain the expected statistical accuracy of the (j>3 measurement as a 
function of B and Dcp statistics. 
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\0 ■ 1- INTRODUCTION 

; 

. Determinations of the Cabbibo-Kobayashi-Maskawa (CKM) [1] matrix elements provide important checks on the 
consistency of the Standard Model and ways to search for new physics. One of the ways to check the consistency of 
^ , the CKM model is measurement of all angles of the unitarity triangle in decays of B mesons. One angle, 0i, has been 

1 ^ ' measured with high precision at modern B factories by the experiments BaBar [2] and Belle [3]. The measurement 
— , of the angle (/)2 is more difficult due to theoretical uncertainties in calculation of the penguin diagram contribution. 

• Precise determination of the third angle, 03, is possible, e.g. , in the decays — > DK^. Although it requires a lot 
' more data than for the other angles, it is theoretically clean due to the absence of loop contributions. Finding the 

^vq 1 way to use the available experimental data most efhciently is therefore essential for the 03 measurement. 
^ ' Various approaches have been proposed to measure 03 in B^ DK^ decays [4-7] by utilizing the interference 

• between b cus and h ucs transitions. The 03 can be measured by studying the interference of B^ D^K^ and 
' B^ — * D'^K^ , which implies the use of neutral D final states common to and Z?". This is provided by using CP 
, eigenstates of D meson (GLW method [4]) or Cabibbo-favored and doubly Cabibbo-suppressed modes (ADS method 

[7]). 

The Dalitz plot analysis of the three-body 13° decay from B^ DK^ processes [8-15] provides today the best 

• measurement of the angle 03. It was proposed by Giri et al. [8] and independently by the Belle collaboration [9]. 

' While Giri et al. mainly discussed a model-independent way of determining 03 using a binned Dalitz analysis of 
'"q^, A'gTT+TT" decay from B^ DK^, the Belle collaboration proposed an unbinned Dalitz plot fit of this decay 

I using the D° AT^TT+TT- model based on a coherent sum of quasi- two-body amplitudes. The latter technique 
has been used so far by BaBar and Belle in their 03 measurement. The most recent BaBar result is 03 = 67° ± 
^ ' 28°(stat) ± 13°(syst) ± ll°(model) using 211 fb^^ data sample [15], the result of the Belle measurement is 03 = 
53° ;j^°(stat)±3°(syst)±9°(model) with 357 fb^^ data sample [13]. The use of Gabibbo-favored mode -> AT^tt+tt^ 
allows to increase the statistics compared to GLW or ADS methods. However, this approach suffers from the model 
uncertainty introduced by the unknown complex phase of decay. Gurrently, this uncertainty is estimated to be 
\J( ' of order of 10°, which is lower than the statistical error. In the future, as the statistics of B decays collected at 
5^ , S-factories increases, this uncertainty will start to dominate the 03 measurement. While the estimations of the model 
uncertainty can also improve, they can only be based on arbitrary variations of the decay amplitude description 
within phenomenological models {e.g. changing the list of resonances or using K-matrix formalism). However these 
models are only approximations of the 3-body dynamics and thus any model uncertainty estimates cannot be treated 
as absolutely reliable. 

It has been shown [8] that the D'^ model uncertainty can be eliminated using the sample of neutral D mesons in 
CP eigenstates. This sample can be obtained at charm facilities, like currently running GLEO-c [16], in the decays 
of -0(3770) into two neutral D mesons. In this paper, we report the results of a feasibility study of such an approach 
using toy Monte-Carlo (MC) simulation based on the realistic description of D decay amplitude. Our main goal is 
the determination of the amount of Dcp decays needed for the model-independent measurement. 

The sensitivity to the angle 03 comes from the interference of two amplitudes producing opposite flavors of neutral 
D meson. Such a mixed state of neutral D will be called D. For example, in the case of i?+ decay, the mixed state 
is D+ = + tbe^^^D^ ■ Here tb is the ratio of the suppressed and favored amplitudes, the total phase between 
the D° and D° is 6*+ = 03 -f (5b, where 5b is the strong phase difference between suppressed and favored B decays. 
Analogously, for the decay of B~ , one can write D- = + rse'^^D" with 9- = —03 + Sb- The ratio tb is 
expected to be of order of 0.1-0.2 [17]. The BaBar (rs = 0.12 ± 0.08 ± 0.03(syst) ± 0.04(model), [14]) and Belle 
{tb = 0.16 ± 0.05 ± O.Ol(syst) ± 0.05(model), [13]) measurements of the tb in Dalitz analysis of from B^ DK^ 
decay confirm these estimates. 
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The Dalitz plot density of D gives immediate information about tb and 9±, once the ampHtude of the D° decay is 
known. The amphtude of the D± decay as a function of DaUtz plot variables = '^/f07r+ ^^'^ ~ ^^k°7t- 

fB± = fD{m^±,m'^) + rBe"^^ fD{m'^,m^±), (1) 

where foifrv^^m"^ ) is an amplitude of the A'^tt+tt^ decay. Here, we neglect possible effects of charm mixing 

and CP violation in D decay. It has been shown [18], that the mixing effects, which arc suppressed in the Standard 
Model, lead to a bias in 03 only at a second order in the small parameters Am£)/r£) and ISTb/^d (where Am^j and 
Ar^) are the mass and decay width differences between the two D mass eigcnstates, F/) is the average decay width 
of mass eigcnstates), and therefore can be safely neglected. 

The — > K^TT^n~ decay model can be determined from a large sample of flavor-tagged K^tt^tt~ decays 

produced in continuum e^e^ annihilation. Once that is known, a simultaneous fit of and data allows to 
separate the contributions of r^, 03 and 5b- However, flavor-tagged D'^ — > i^^Tr+Tr" decays only provide information 
about the absolute value of the amplitude. The phase term can only be obtained based on some model assumptions 
(such as Breit-Wigner amplitude dependence of the underlying resonances) which results in the model uncertainty of 
the 03 measurement. 

In the original publication by Giri et al. [8], an approach based on the binned Dalitz plot analysis was proposed 
which allows to extract 03 without model uncertainties. The idea of the approach is based on the fact that, to measure 
03, complete knowledge of the phase term of Jd amplitude is not necessary. In the D'^ —> Kgir'^Tr' decay amplitude, 
we emphasize the absolute value |/| and phase S^: 

fD{ml,ml) = 1/1,(774, m2_)|e'^^-("^"-) (2) 

Since the phase of fs is arbitrary, it can be expressed as 

/b± - 1/1,(7771, 777^)1 +rse±'^^^+'^- 1/^,(7772^, 777^^)1 exp(i5p(7n|, 777^) -75^(^2^,7772^)), (3) 

i.e. only the difference of phase terms SDirn^jiri^) — (J/, (777 2^, 777^) between symmetric points of the phase space is 
relevant. 

This difference can be obtained either by studying the decays of neutral D in CP eigenstate or treating them as free 
parameters in the extraction of the 03 from the B decays (see Section 5). Decays of D mesons in CP eigenstate to 
KgTr~^TT^ can be obtained at e+e^ machine in the process, 6.17. e+e~ — *■ ?/'(3770) — > DD, where the other (tag-side) D 
meson is reconstructed in CP eigenstate, such as K~^K~ or KgOj. The binned Dalitz plot analysis allows to determine 
03 without the need to parametrize the dynamics of the D*^ —> Kc;Tt~^tt~ decay {e.g. by introducing resonances) and 
thus allows for a truly model independent measurement. 

Different approaches to the extraction of Sjj phase and different Dalitz plot binnings arc tested in our feasibility 
study and compared to the model-dependend approach performed in an unbinned fashion. We consider only the 
approaches described in [8]. However, another possibility exists to access the phase information, which involves the 
process e+e^ ^{3770) DD with both neutral D mesons decaying to KgTr~^Tr~ final state [19]. This approach 
requires to study the correlations of the two Dalitz plots. 

In the next section, we follow the technique proposed in [8] , to introduce the notation which will be used in the rest 
of the paper. 



2. BASIC IDEA OF THE MODEL-INDEPENDENT TECHNIQUE 

The amplitude of A'^tt+tt^ decay from CP eigenstate of neutral D is 

fcP±{mX,m'i) = ■-i=[/D(n7+,?77?.) ± /d(777?., 777+)], (4) 

for CP-even and CP-odd states of neutral D, respectively. 

Consider the D° — > A'^tt+tt^ Dalitz plot is divided into 2JV bins symmetrically over exchange 777^ m^. The 
bins will be denoted by the index i ranging from —Af to Af (excluding 0); the exchange <-> rril_ corresponds to 
the exchange i ^ — i. The number of events in the 7-th bin of the Dalitz plot, given by the region Pj, is 

K, = Ad[ \fDiml,ml)\^dmldml = AdF„ (5) 

JV, 

where A/, is the normalization that depends on the luminosity integral and definition of amplitude fo- 
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Similarly, as follows from Eq. 1, for i^^Tr+TT decays from DK^ mode the number of events in z-th 

bin is 



/ \fD{rn\,m^zf.)\'^dm^_dm\+r% I \fD{ml^,m\)\^dm^_drn\ 



+2rBi?e (e*''± / fD{m\,m^_)f}j{rn^_^,in\)dm^_dm 



(6) 



= Ab [F±, + 2rBRe{e''>±Gi) + r|F^,] 

= Ab [F±i + 2rB{cos9± Re d ~ sin6'± Im Gi) + r%F^i] . 

Thus in order to extract information about 9±, one needs the values of F±i, which are obtained from the decays of 
flavor tagged D mesons, and also the values of Re Gi and Im Gi. 
The number of events in z-th bin for CP eigenstate of neutral D is 



M^=AcP±j \fcp±{ml,ml)\''dmldml 
= 2AcP±{F,±2Re G., + F.,). 



(7) 



Substituting the integrals Fi from Eq. 5 and introducing the parameters Ci and Si which depend on the phase term 
4>{m'^ , m?_ ) , we finally obtain the following system of equations which relates the number of events in bins of flavor 
D, Dcp and D from B decays: 



hB[K±i + t'^bK^i + 2rB^/K~K~{ct cos(±(?!)3 + 5b) - Si sm{±(j)3 + Sb))], 



(8) 



M± = hcP± [K, + ± 2^K~KZ,Ci\ , (9) 

where hp and hcp± are the normalization factors which take into account different branching ratios, luminosity 
integrals and detection efficiencies for B, and Dcp± decays, 

\fD{m\,mt)\\fD{mi,m'l_)\ cos(5£)(to^, m^) - 5D{rri^-,m\))dm^_dm\, (10) 



and 



|/n(m|,mi)||/£,(TOi,77i^)|sin((5n(m^,mi) - 5D{mt,m\))dmtdm]^. (11) 



Note that Ci = c-i and Si = —S-i. The meaning of these parameters is sine and cosine of the phase difference, 
averaged over the bin area with the weight proportional to m?_)||/D(m?_, m^)|. It is clear that sf + cf < 1. 

The parameters Ci and Si are in general unknown, they can be determined either from Dcp decays, or by leaving 
them as free parameters in the fit to D Dalitz plot from B^ DK^ decay. Depending on the way these parameters 
are obtained, several approaches of 03 extraction are possible. Below we discuss these approaches, together with the 
toy MC simulation results for each of them. 

In our toy MC studies, we assume that the statistics of flavor-tagged D mesons is practically infinite, therefore, the 
absolute value of the fn amplitude is known to high precision. This is provided by the large cross-section of D*^ 
meson production in e+e~ — > cc process. We also assume that the number of Dcp decays is much larger than the 
number of B decays, due to the fact that the cross-section of e+e~ — > il;{3770) D'^D^p production is approximately 
two orders of magnitude larger than that of e+e~ T(45) BBdk process. 



3. APPROACH WITH TAGGED Dcp AND CONSTRAINT ON Si 



It is clear that for the most precise measurement of (f)^ given the limited statistics of B decays, all the available 
information has to be used. In particular, the unknown quantities Ci and Si have to be extracted from the higher- 
statistics Dcp data. 
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The coefScient Ci can be obtained directly from the and Dcp decays. For example, if the numbers of Dcp+ and 
Dcp- decays are proportional to the branching fractions of ?/'(3770) — > DDcp+ and -0(3770) — > DDqp-, respectively, 
i.e. if hcp+ = hcp-, 

c.4^^-^^^ + ^- (12) 

(see Eq. 9). However, c,; can be extracted even in D decays of only one CP parity. 

If the binning is fine enough so that both f{m\, m"L) and f(m?L, m^) are approximately constant over (m^, m^) G 
T>i, the expressions for a (10) and Si (11) reduce to Ci = cos {Spi{m'^,m'^) ~ (5_D(m?_, m^)), Si = sin((5£)(m^, to?.) — 
5_D(m?-, TO^)). In such a case, Si can be directly obtained as Si ~ ±-\/l — c^. The sign of Si can be chosen either 
based on the continuity requirement (it should vary smoothly for adjacent bins) or under the assumption that the D° 
decay proceeds via two-body amplitudes. For simplicity, in our study we take the "correct" sign (calculated from 
the input Z?" model). 

Once the coefficients and Si are determined, the parameters r^, 03 and 8b can be obtained from B decays by 
minimizing the negative logarithmic likelihood 



-2iog/: = -2iog/:+ - 2iog/:_, (la) 



where 



-2l0g/:±= — — ' ' -^V"'"-'v^»— v-ro . , -^^nu ^^^^ 



{Nf - hB\K±^ + r\K^^ + 2rB^K^KZ^{c^ cos(±03 + fe) - s^ sin(±03 + fc))])^ 



is the negative logarithmic likelihood for one B flavor (5+ or i?^), N'^ is the number of events in z-th bin for 
corresponding B flavor. Here, Gaussian behavior of the number of events in each bin is assumed. 

To avoid possible biases due to non-constant /(to^,to?_), the bin size should be relatively small. The size of the 
bin is determined by the width of the most narrow structure in the Dalitz plot; in our case, it is the lo meson with 
F = 8.4 MeV/c^; therefore, typical scale of phase variation is Ato^ = 0.013 GeV^/c*. However, since we assume that 
the number of B decays is much smaller than the number of Dqp decays, the number of B events per bin for such 
a fine binning can be small. This can affect the fit procedure since we expect Gaussian behavior of the number of 
events. Therefore, for our study we take relatively rough binning for the fit of B decays, and introduce a sub-binning 
of each of "i?-bins" into several "Z?-bins". As follows from definitions (10), (11) and (5), the parameters Ci and si for 
" _B-bins" can be calculated from finer " D-bins" as 



Cj,Q yjKi^aK-i-a/ ^jKiK^i, (15) 



Si = 'Y^Si,a^ Ki,aK^i^^al \/ KjK^i, (16) 

where the index a corresponds to the sub-binning of the bin TDi. 

A toy MC study has been performed to obtain the accuracy of 03 measurement as a function of B and Dc p decay 
statistics, and to investigate the influence of the finite bin size on the measurement result. The Dalitz plot binning 
used for our MC study is shown in Fig. 1. The phase space Am^ = rrii^ — rr?_ e (—3 GeV^/c^,3 GeV^/c^) and 
TO^^ = M|, — — M\ — — vr?_ £ (0, 2 GeV^/c^) was uniformly divided into rectangular bins, symmetrically 
over the exchange of the sign of the final state pions. The bin size should certainly affect the statistical precision of 
the measurement. Since the structure of the D° KgTT~^Tr~ decay matrix element is quite complex, the variations 
of the event density due to CP violation change rapidly across the phase space, and, as the bin size becomes larger, 
the average effect of CP violation on each bin is reduced. Three different binnings were studied; depending on the 
number of bins in half of the phase space they are referred to as 3x3 (the corresponding number of bins which enter 
the half of the kinematically allowed phase space is Af = 8), 5x5 {Af = 19) and 7x7 {Af = 32) binnings. 

The parameters rs and 5b were taken to be 0.2 and 180°, respectively; 03 was scanned in a range from 0° to 180°. 
For the D° decay model, the recent Belle measurement was used [12]. The efficiency variations over the phase space 
and background contribution were not taken into account in the simulation. The statistical precision of the method 
certainly depends on the actual D'^ model, therefore, the real experimental precision can somewhat differ from our 
estimate. The error is also inversely proportional to the magnitude of the opposite D flavor contribution. The 
estimates of this value vary in the range 0.1-0.2, the current measurements of tb [13, 15] agree with these numbers. 
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FIG. 1: Binning of the Z?" — > K^-k^tt Dalitz plot. 3x3 (left), 5x5 (center) and 7x7 binning (right). 5x5 sub-binning of one of 
the bins of the central plot is shown; the sub-binning is used for Dcp Dalitz plot. 

although the error is still large. If the actual value of differs from our expectation, the results of our study should 
be scaled accordingly. 

The statistical error of the (t>3 measurement as a function of (jj^ for the different binnings is shown in Fig. 2. For 
each (j>3 value, 200 pseudo-experiments with the statistics of 10^ — > DK^ decays were generated. The (f>3 error 
was obtained from the spread of the reconstructed values. To keep the contribution of finite Dcp statistics small, we 
used the large number of Dcp decays (2 x 10^ decays of each CP-parity). The 03 accuracy depends only weakly on 
the values of 03 or Sb since both c,; and Si are obtained from the same Dcp statistics and have comparable accuracies. 
There is no significant difference in the statistical accuracy for 5x5 and 7x7 binnings, while in the case of 3x3 binning 
the accuracy is worse by 30-40%. We have used the 5x5 binning for all further studies (excluding the study of the 
approach without Dcp decays, sec Section 5). The precision of 03 extracted from a model-dependent unbinned Dalitz 
plot fit with the same input parameters is also shown for comparison. The unbinned fit is the optimal technique from 
the statistical point of view and represents the best limit of the statistical precision. The resolution in that case is 
approximately 30% better than for the case of binned approach. The optimal choice of binning can possibly improve 
the accuracy of the binned technique. We leave that question for the future studies. 

The results of the toy MC study of 03 bias due to finite D-bin size are shown in Fig. 3. The angle 03 was scanned 
from to 180°, other parameters were taken to be tb ~ 0.2, 5b ~ 180°. Points show the 03 bias for no sub-binning, 
3x3 sub-binning and 5x5 sub-binning (shown in Fig. 1 for one of the bins). For our further study, we have used 5x5 
sub-binning, for which the bias does not exceed 2°. In a real experiment the bias can be corrected, e.g. based on 
the quasi two-body model. Since the correction is comparatively small, the uncertainty introduced by a model-based 
correction should be negligible. 

The results of the study of 03 sensitivity dependence on the number of reconstructed DK^ and Dcp decays 

are shown in Fig. 4. The graph shows the 03 error for tb = 0.2, 03 = 70°, Sb = 180°. For different values of r^, 
the error scales inversely proportional to tb- The weak dependence of the statistical error on the phases 03 and 6b 
should also be noted, therefore, the actual 03 error can differ by 30-50% depending on the true values of these phases. 
The dependences on both B and Dcp statistics show a square root scaling. For comparison, the statistical error of 
03 extracted by model-dependent unbinned Dalitz plot fit with the same input parameters is shown with the dotted 
line. 

According to our estimations, the proposed super-B factory with its design integrated luminosity of 50 ab^^, would 
allow a measurement of 03 with accuracy below 2° . To keep the contribution to the statistical error from the finite 
sample of Dcp decays below that level, one needs about 10** tagged Dqp decays, corresponding to approximately 
10 fb~^ of integrated luminosity at V'(3770) peak. In the following sections we consider possible alternative approaches 
of the binned analysis and compare their precision with the approach described above, which we have found to be the 
most optimal. 

4. APPROACH WITHOUT Si CONSTRAINT 

Constraining Si from Ci introduces a slight model dependence to the measurement result, since the bin size becomes 
dependent on the width of the Dalitz plot structures of D° decay. In principle, the system of equations is solvable 
even without this constraint. In this case, there are SAf equations (8, 9) for 2J\f + 6 unknowns (q, Si, tb, 03, 5b, 
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FIG. 2: 03 error as a function of (jta, for different Dalitz plot 
binnings. Results of the toy MC study with 2 x 10* — > 
DK^ decays, 4 x 10^ Dcp decays, rs = 0.2, = 180°. 



FIG. 3: 03 bias due to finite size of the bin as a function of 
03, for the cases of no sub-binning (circles), 3x3 sub-binning 
(triangles), 5x5 sub-binning (squares). 




FIG. 4: Statistical error of 03 as a function of number of reconstructed B — » DK decays (open circles) and Dcp decays 
(dots). Results of toy MC study with rp = 0.2, 03 = 70°, 5b ~ 180°. 03 error from model-dependent unbinned Dalitz plot fit 
with the same input parameters is also shown (dotted line). 
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FIG. 5: <j)3 error as a function of 03 for relaxed Si (triangles); 
and Si constrained from Ci (squares). Results of the toy MC 
study with 2 x 10"* DK^ decays, 4 x 10^ Dcp decays, 

tb ~ 0.2, Sb ~ 180°, 5x5 Dalitz plot binning. 



FIG. 6: 03 error as a function of 03 for the model- 
independent approach without Dcp decays (triangles) and 
with Dcp decays and Si constraint (squares) . Results of the 
toy MC study with 2 x 10^ DK^ decays, tb = 0.2, 

6b = 180°, 3x3 Dalitz plot binning. 



hs, hcp±). CoefBcients Ci are still directly obtained from the Dcp decays, while Si are mainly constrained by the B 
decays. Therefore, we expect the accuracy of the cos(±(/)3-|-i5b) determination to be higher than that for sm{±(j)^ + 6B)- 
In other words, the accuracy of (j)^ determination should strongly depend on the values of (p^ and 5b- The 03 error 
as a function of (j>z for rs = 0.2 and 5b = 180° is shown in Fig. 5 (triangles). This plot was obtained by toy MC 
procedure with 10^ DK^ decays of each flavor, and 2 x 10^ Dcp decays of each CP-parity. For comparison, 

the 03 error is shown obtained with the approach where Si is fixed from c;. As follows from Fig. 5, the precision of the 
approach under study is comparable with the optimal strategy precision only in the narrow range of 03 values. The 
position of this region depends on the value of the strong phase 5b, therefore, the precision of the 03 determination 
in this approach can not be guaranteed. 



5. APPROACH WITHOUT TAGGED Dcp 



Even if the Dcp decays are not available, it is in principle possible to extract 03 in a model- independent way using 
only DK^ and flavor-tagged D decays. For Af pairs of bins, there are 4J\f equations (8) {2JV for each flavor of 

B), and 2Af + 4: unknowns (q, Si, tb, 4>3, 5b, fiB)- Therefore, the system is solvable for M > 2. Technically (if the 
number of equations is greater than the number of unknowns), it can be done by minimizing the negative logarithmic 
likelihood (Eq. 13, 14) with Ci and coefficients taken as free parameters. 

Practically, however, this approach offers too low sensitivity for any reasonable B decay statistics. A toy MC study 
has been performed with 10^ B^ — > DK^ decays of each fiavor. The parameters and 5b were taken to be 0.2 and 
180°, respectively; 03 was scanned in a range from 30° to 140°. The results of the toy MC study are presented in 
Fig. 6 for 3x3 binning. Statistical precision with Dcp data and Si constraint is also shown for comparison. It can be 
seen from the plots, that the accuracy of the 03 determination is quite poor compared to the method which makes use 
of Dcp decays. Reaching the accuracy of 2° would require about 5 x 10^ reconstructed B^ DK^ decays, which 
corresponds to the currently unrealistic luminosity integral of 5000 ab~^ at B-factory. 

Taking finer 5x5 binning could slightly improve the sensitivity, however, due to large number of free parameters, 
the fit converges poorly even with a number of B decays as large as 10^. Variations of the initial parameters in the fit 
can lead to significant changes of the fit result. Binning finer than 5x5 has not been investigated due to an internal 
limitation on the maximum number of free parameters in the software used for minimization. 
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6. CONCLUSION 

We have performed a toy Monte-Carlo study of a model-independent way to measure the angle 4>3 of the unitarity 
triangle, suggested by Giri et al. [8]. The method involves — > DK^ decays with neutral D decaying to i^^Tr+Tr" 
final state, together with the decays of D mesons in CP eigenstate to the same final state. Different approaches to the 
extraction of (p^ are considered. Although technically the model-independent measurement of <^3 can be performed 
without decays of D in CP eigenstate, the statistics of — > DK^ decays required to obtain a reasonable 03 accuracy 
in this approach is unrealistically high. 

We have shown that for the optimal strategy (which involves constraining Si coefficients from Ci), the accuracy 
of the model-independent binned approach is only 30% worse compared to the unbinned technique. Possibly the 
choice of the optimal binning {e.g. using non-rectangular bins) can improve the accuracy of the model-independent 
measurement. The statistical error of 03 does not depend significantly on the 03 and strong phase values. 

According to our estimations, the proposed super- i? factory [20] with its design integrated luminosity of 50 ab~^, 
would allow a measurement of 03 with statistical accuracy below 2° . To keep the contribution to the statistical error 
from the finite sample of Dcp decays below that level, around 10^ tagged Dcp decays are needed, corresponding to 
approximately 10 fb"-'^ of integrated luminosity collected at -0(3770) peak. This task can be accomplished by CLEO-c 
experiment at CESR collider [16] and by upgraded BESIII/BEPCII complex at Beijing with design luminosity of 10'^'^ 
cm-2s-i [21]. 

The systematic errors due to variations of the detection efficiency or background density over the Dalitz plot 
and other related factors were not considered in this study. They will certainly depend on the specific detector 
parameters, but the current experience with model-dependent Dalitz analyses at B factories allows us to believe that 
these contributions will also decrease with increasing statistics. 
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